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ABSTRACT 



We study the effect of local stellar radiation and ultraviolet background radiation 
(UVB) on the physical properties of damped Lyman-a systems (DLAs) and Lyman 
■ limit systems (LLSs) at z = 3 using cosmological SPH simulations. We post-process 

, our simulations with the Authentic Ray Tracing (ART) code for radiative transfer of 

local stellar radiation and UVB. We find that the DLA and LLS cross sections are 
significantly reduced by the UVB, whereas the local stellar radiation does not affect 
them very much except in the low-mass halos. This is because the clumpy high-density 
O ■ clouds near young star clusters effectively absorb most of the ionizing photons from 

young stars. We also find that the UVB model with a simple density threshold for 
C/3 . self-shielding effect can reproduce the observed column density distribution function 

^ ■ of DLAs and LLSs very well, and we validate this model by direct radiative transfer 

calculations of stellar radiation and UVB with high angular resolution. We show that, 
^SJ ■ with a self-shielding treatment, the DLAs have an extended distribution around star- 

^ I forming regions typically on 10 — 30 kpc scales, and LLSs are surrounding DLAs 

on ~ 30 — 60 kpc scales. The DLA gas is less extended than the virial radius of 
0^ ■ the halo, and LLSs are distributed over the similar scale to the virial radius of the 

host halo. Our simulations suggest that the median properties of DLA host haloes are: 
^ I Mh = 2.4 X 10^° Mq, SFR= 0.3 Mgyr"!, = 2.4 x 10^ Mq, and Z/Zq ^ 0.1. About 

PsJ ■ 30 per cent of DLAs are hosted by haloes having SFR =1 — 20 Mq yr~^, which is the 

typical SFR range for LBGs. More than half of DLAs are hosted by the LBGs that 
are fainter than the current observational limit. Our results suggest that fractional 
contribution to LLSs from lower mass haloes is greater than for DLAs. Therefore the 
median values of LLS host haloes are somewhat lower with Afh — 9.6 x 10^ Mq, 
SFR= 0.06 Mgyr-i, 6.5 x 10^ Mq and Z/Zq = 0.08. About 80 per cent of 

total LLS cross section are hosted by haloes with SFR < 1 Mq yr~^, hence most LLSs 
M ■ are associated with low-mass halos with faint LBGs below the current detection limit. 

Key words: radiative transfer - galaxies: evolution - galaxies: formation - galaxies: 
high-redshift - galaxies: ISM - methods: numerical 



1 INTRODUCTION 

Interplay between neutral hydrogen gas and star for- 
mation in galaxies is a crucial issue for understand- 
ing galaxy formation and evolution. Absorption systems 
in the spectrum of QSOs are good tracers of neu- 
tral hydrogen, and large surveys such as the Sloan 
Digital Sky Survey (SDSS) provide abundant data on 
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Hi g a s at high redshift (IProchaska. Herbert-Fort fc Wolfel 



l2005l : IProchaska fc Wolfel 120091 ') . The absorption systems 
with A^Hi ^ 2 X 10'^° cm~^ show the damping wing 
in the absorption profile, and are c alled Damped Lya 
systems (DLAs) (|Wolfe et all Il986l ). In addition, re- 
cently the D LAs has been studied even by GRB af- 
tergrow (e.g.. [ Prochaska et all l2007al : iLedoux et al. I l2009l : 
iRau et all l20ld ). It is known that DLAs dominate the 
ne utral gas conten t of Universe between z ~ — 
5 (| Wolfe et all 1 19861 : IStorrie-Lombardi. McMahon fc IrwinI 
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Il996l : IStorrie-Lombardi fc Wolfd [ioool ') . The global neutral 
gas content measured from DLAs declines with cosmic 
epoch in concert with the grow th in cosmic stellar mass 
(IStorrie-Lombardi fc Wolfd I2OOOI ). Therefore, it is generally 
thought that DLAs are strongly correlated with cosmic star 
formation. However, in spite of numerous observational sam- 
ple of DLAs, a deep understanding of DLAs and host galax- 
ies remains elusive. 



by recent progress in observations (e.g., Wolfe et al. 


1995 


: IStorric-Lombardi & Wolft 


2OOOI: Rao & Turnshek 


200c 


: Peroux et al. 20031: 


Chen & Lanzettal 


2002; 


Prochaska et al.1 20051: Prochaska & Wolfd 


200c ; 


Prochaska, O'Meara & Worseck 201G|: Wolfe et all 


2008; 


Noterdaeme et al. 12009*) and simulations (e.g., Katz et al. 


199e 


: Gardner et al., 1997a,b: Hachnolt, Steinmetz & Ranch 


1998 


: iMaller et al.1 I2001I: iGardner et al.1 I2OOII: 


Razoumov et al.l l2006al: 


Nagamine. Springel & Hernauist 


2004b a; Nagamine et al 


20071: iPontzen et al.l I2OO8I: 


Barnes & Haehneld l2009l: 


Tescari et al.1 l2009l: ICenI l2010h. 



Thanks to these works, our understanding on the relation- 
ship between DLAs and host galaxies have improved. For 
example, it has been suggested that a significant fraction of 
DLAs are distributed in the same dark matter halos as the 



Lyman break galaxies (LBGs; e.g., [A dclbcrgcr ct al. 2005' 



ICooke et al]|2006l : iRafelski. Wolfe fcChcn.,2011, ). Theoret- 
ically, [Ce3 (|2010D examined the number fraction of DLAs 
related t o LBG s using cosmological AMR simulations. 
iLee et al] (|201ll ') estimated the DLA-LBG cross-correlation 
function using cosmological SPH simulations, and showed 
th at their resul t agree s well with the observational estimate 
of ICooke et al] (|2006l ). 

A connection betwe en DLAs and Lvq em itters (LAEs) 
has also been suggested. iRahmani et aL (|2010l ) investigated 
Lya emission from 341 DLAs at 2 ~ 2.86, and concluded 
that the overall DLA population seems to origi nate from the 
low lu minosity end of LAEs at high redshift. iBarnes et al.l 
(|201ll ) simulated three dimensional Lya radiative transfer 
(RT) for high-redshift galaxies in cosmological SPH simula- 
tions. They showed that DLA host galaxies show extended 
Lya emission around high-column density region. The Lya 
emission is extended over several arc seconds, and has the 
spectral width of several hundred km s~^. A combination 
of these observational and theoretical works seem to suggest 
that DLAs are associated with a wide variety of sources, 
such as LBGs and LAEs. 

For a better understanding of DLAs, it is quite impor- 
tant to investigate various feedback effects in galaxies which 
can drastically change the state of gas. Multiple supernovae 
can cause galactic wind, and hence change gas distribution 
largely in galaxies. iNagamine et alj (l2004bf ) investigated the 
effect of galactic wind on the H I column density distribution 
function /(A^hi), and found that strong galactic winds with 
a speed of ~500kms~^ is n eeded to reproduce t he observed 
/(Afni) at high A^hi values. iTescari et all (|2009t ) also found 
that strong galactic winds with speed of ~600kms^^ are 
needed in order to reproduce the observed f{Nm) and the 
evolution of Hi mass density with redshift. Several authors 
have suggested that galactic outflows can produ ce the large 
veloc i ty dispersion as observed i n DLAs (e.g., iRazouniovl 
I2OO9I : iHong et al.ll2010l : FCenlbOlOt ). 

Gas inflows from surrounding IGM can also have 



strong impact on the absorption systems in high-2 galax- 
ies. Recent simulations have revealed a new theoreti- 
cal paradigm for galaxy evolution that a large amount 
of gas penetrate deep in side dark matte r haloes as 
cold, filame n tary streams ( Katz et al.l I2OO3I; ores ct al.l 
20051 . I2OO9I: iBiimboim fc Dek cl 2003; Dckcl & Birnboim 



2ooe 



Ocvirk. Pichon fc Tevss icr 2008; Brooks ct al. 20o3; 



IDekel et al.ll2009l l. iFaucher-Gigucrc fc Korea (,2011. ) studied 
the neutral hydrogen content in the cold accretion to LBGs 
using SPH simulations with UVB RT. They showed that a 
part of gas could be DLAs , however the covering fraction 
was quite small. Similarly, iFumagalli et al.l (|201ll ) showed 
the neutral hydrogen gas content in cold stream by cosmo- 
logical AMR simulations and RT calculation of local stel- 
lar source and UVB. They find that most of the gas in 
cold streams is ionized, however it can survive as Ly man- 
limit systems (LLSs), which is defined by column densities 
of 10" < A^Hi < 10^°-^ atoms cm'^ 

The observations of DLA metallicity can also give im- 
portant constraints on feedback processes, as the interstel- 
lar ga s can be enriched by super novae in the early Uni- 
verse. iMoUer. Fvnbo fc Fall (|2004l ) presented the metallic- 
ity distribution of DLAs and suggested that the DLAs 
with Lya emis sion tend to have higher metallicities of 
Z > 10"i-^ Zf^^. lLedoux et al.1 (|2006l ') presented the velocity 
width Au-metallicity relation of DLAs from the line width 
of low- ionization metals, and showed that the metallicity is 
linearly proportional to Av with a large dispersion. They 
also suggested that the slope of fitting function for Aii- 
metallicity relation is the same for the sample at z > 2.43 
and 2 ^ 2.43, but the medi an metallicity decreases with 
increasing redshift. Moreover. iProchaska et al] (|2007lJ ) pre- 
sented the metallicities of a lar ger DLA samp l e with a large 
dispersion similar to that of iLedoux et al.l (|2006l ). Theo- 
retically, recent numerical simulations which include metal 
enrichment by supernovae have succe ssfully produced the 
observed metal distribution o f DLAs (|Pontzen et al.|[2008l : 
iTescari et al]l2009l : [Cenll2010l '). 

Another important feedback on DLAs is radiation 
from UVB and local stellar sources. The UV radia- 
tion can ionize hydrogen gas, and largely reduce the 
neutral fraction of ISM. Ea rlier Katz et all (| 19961 ') and 
iHaehnelt. Steinmetz fc RauchI (|l998l ) included the self- 
shielding correction for UVB radiation to reproduce DLA 
statistics in there numerical simulations. Thereafter, recent 
simulations have reproduced various observation data 
of DLAs reasonably wel l by c o nside ri ng the effect of 
UV R T (Razoumov et al.1 l2006al. I2OO8I: iKohler fc Gnedinl 
20071 : iPontzen et al.1 I2OO8I: iNagamine. Choi fc Yaiimal 



20ld: lAltavet al.1 I2OI1I: iFaucher-Giguere fc Kere^ |201~ 



Fumagalli et al.l I2OIII : iMcQuinn. Oh fc Faucher-Giguere 1 
20111 ). Most of them have focused on the RT of UVB, and 
found that the self-shielding effect is crucial for m a tching 
the observed f{Nm) of DLAs. INagamine et all (|20ld ) 
showed that the simple self-shielding model without RT 
calculations could reproduce the observed /(A^hi) very 
well e xcept for the very high A'hi values. lAltav et al.l 
(|201ll ) also reproduced the observed /(A'hi) in the range 
of lO^'^ < Ahi < 10'^'^ atoms cm~^ by post-processing 
cosmological SPH simulations for the UVB RT. 

On the other hand, the fiux of stellar radiation near the 
star cluster can be much higher than the local UVB, hence 
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stellar RT can be important for DLAs and LLSs. However, 
it has been ignored in many simulations. Some authors have 
suggested that the UV radiation from loc al stellar sources 
can be important for absorption systems (|Miralda-Escudel 
I2OO5I : ISchavll2006l ). however, there have not been exten- 
sive studies of direct RT cal culations that examine d the ef- 
fect of local stellar sources. iNagamine et al.1 (|2010l ') showed 
that the f{Niii) at z = 3 does not change very much, 
based on a direc t RT c alculation from young stellar sources. 
iFumagalli et al.l (|201ll ) carried out the RT calculations of 
local stellar sources and UVB for seven high-resolution sim- 
ulated galaxies, and found that the local stellar radiation 
changed f{Nm) at Nm ~ 10^® - 10^^ cm"^ by ~0.5 dex 
to the other model without local stellar radiation. Their 
simulation result is in good agreement with observation at 
logAi'm > 21, however it underestimates f{Nm) at lower 
A^Hi range of LLS. 

The stellar radiation flux onto a gas cloud varies de- 
pending on the distance between gas and stars and the 
type of SED. Therefore in a more realistic inhomogeneous 
medium, precise RT calculations are needed to study the 
effects of stellar radiation on absorption systems. Further- 
more, since the distribu tion of stars and gas varies among 
different galaxies (e.g., lYaiima et al.l l2011al ). we need to 
carry out the RT calculation for many galaxies to dis- 
cuss the statistical data such as f{Nm)- In the present 
work, we perform precise RT calculations with ~ 10* an- 
gular bins for each star particle in cosmological SPH sim- 
ulations, and examine the distribution of neutral hydrogen 
gas for a few hundred galaxies. Our RT code is based on 
the ray-tracing method, which can estimate the ionization 
stru cture more accurately than the Monte Carlo met hod 
that IFumagalli et al.l (|201ll ) used (see llliev et aL I l2006l . for 
a comparison of RT methods). In addition, the number of 
galaxy sample processed with RT in our wo r k is m uch larger 
(A^ > 100) than that of iFumagalh et al. I (|201ll . N ~ 8). 
Thus, our calculation allows us to investigate the effect of 
RT more reliably on the statistical quantities of DLAs and 
LLSs, such as /(Ahi) and the mass dependence of cross sec- 
tions. Moreover, by direct RT calculations o f UVB, we verify 
the se lf-shielding model of UVB derived bv lNagamine et al] 
(|2010l ). We also discuss the physical quantities of host galax- 
ies of DLAs and LLSs, such as the halo mass, star formation 
rate, stellar mass and metallicity. 

Our paper is organized as follows. We describe our sim- 
ulations and our approach for RT calculations in Section [51 
In Section [3l we present our results, and show the effect of 
stellar and UVB radiation on cross sections of DLAs and 
LLSs, and on f{Nui)- In Section [3.51 the cumulative abun- 
dance of DLAs is shown. In addition, we show the typical 
physical quantities of DLAs and LLSs host galaxies. Finally, 
we summarize our main conclusions in Section [3| We focus 
on redshift 2 = 3 in this paper. 



2 MODEL AND METHOD 
2.1 Simulations 

We use an updated and modified version of the Tree- 
particle-mesh (TreePM) smoothed particle hydro dynamics 
(SPH) code GADGET-3 (originally described in ISpringell 



iooi). The SPH calculation is performed based on the en - 
tropy conservative formulation ([Springel fc Hernquistl200^ ') . 
Our fiducial code includes radiativ e cooling by H, He, 
and metals (|Choi fc Nagaminll2009l '). star formation, su- 
pernova feedback, a phenomenological model for galac- 
tic winds and a sub-resolu tion model of multiphase ISM 
([Springel fc Hernauistl2003l l . We also include the heating by 
a uniform UVB, which we will discuss more in Section [2. 2. 2 1 

In this multiphase ISM model, high-density ISM is 
pictured to be a two-phase fluid consisting of cold clouds 
in pressure equilibrium with a hot ambient phase. Cold 
clouds grow by radiative cooling out of the hot medium, 
and this material forms the reservoir of baryons available 
for star formation. The star formation rate (SFR) is esti- 
mated for each gas particle that have densities above the 
threshold density, and the star particles are spawned statis- 
tically based on the SFR. For the s tar formation model, the 
"Pressure SF model" is being u sed ([Schave fc Dalla Vecchial 
[2OO8I : [Choi fc Nagamind [201oh . This model estimates the 
SFR based on the local gas pressure rather than the gas 
density, and prevents artificial fragmentation of the gas. 

The simulations used in this paper also uses the 
"Multicomponent Variable Wind" model developed by 
[Choi fc Nagamind ([201lh . This wind model is based on both 
the energy-driven wind and the mo mentum-driven wind 
[Murrav. Ouataert fc ThompsorJ ([2005h . and the wind speed 
depends on the galaxy stellar mass and SFR. It gives more 
favorable results when compared to the observations of cos- 
mic C IV mass density and IGM temperature than the pre- 
vious model wit h a constant wind speed . To enable this 
new wind model, [Choi fc Nagamind (|201lf ) implemented an 
on-the-fly group-finder into GADGET-3 to compute galaxy 
masses and SFRs while the simulation is running. The 
group-finder, which is a simp lified variant of the SUB- 
FIND algorithm developed by [Springel. Yoshida fc Whitd 
([200lh . identifies the isolated groups of star and gas parti- 
cles (i.e., galaxies) based on the baryonic density field. The 
detailed procedure of t his galaxy grouping is described in 
[Nag aminc et all ([2004dl 'l. 

The simulations used in this paper are performed with 
2 X 144"^ particles. The DM and gas particles masses are 
1.97 X 10^ and 4.06 x 10*^ /i'^Mq, respectively. The comov- 
ing gravitational softening length is e = 2.78h~^ kpc in a 
comoving (10/i^^Mpc)^ calculation box. 

The adopted cosmolog ical parameters are c onsistent 
with the WMAP results ([Komatsu et al.l [201J): Ho = 
72kms-^Mpc~^ {h = 0.72), D.m = 0.26, = 0.74, 
Qb = 0.044, as = 0.80, and Us = 0.96. 



2.2 Radiation Transfer 

2.2.1 Stellar Radiation 

To study the effect of stellar radiation on DLAs and LLSs, 
we compute the stellar radiation transfer and the ioniza- 
tion structure of gas in each dark matter halo by post- 
processing the simulation output. The RT scheme used 
in this paper is the Authentic Radiation Transfer (ART) 
method, and the tr e atment is basically the sa m e as i n 
Yaiima et al] ([2009l'l: [ Yaiima. Oioi fc Nagamind (|2011a[ ): 
Yaiima. Umemura fc Mori ( 2011bh . Here we briefiv summa- 



rize the procedure of RT calculation. 
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First we set up a uniform grid around each dark matter 
halo with a grid cell size equal to the gravitational softening 
length of the simulation, and translate the SPH gas informa- 
tion into a gridded data. As for the scatte ring of photons, w e 
employ the on-the-spot approximation (|Osterbrockl 1 1989f ) . 
in which the scattered photons are assumed to be absorbed 
immediately on the spot. In this work, the RT equation is 
solved along Tig rays (~10*~^) with a uniform angular res- 
olution from each source, where ng is the grid number on 
a side of the calculation box of each galaxy. We shoot the 
radiation rays in a radial fashion from each star particles. 
The number of ionizing photons emitted from the source 
stars is computed based on the theoretical spectral energy 
distribution (S EP) given by the po pulation synthesis code 
PEGASE v2.0 (|Fioc fc Roccalll997l ). We consider the effect 
of age and metallicity of the stellar population by interpo- 
lating the t able gene r ated f rom the result of PEGASE. We 
assume the ISalpeteil (|l955h initial mass function with the 
mass range of 0.1 — 50Mq. 

In addition, we c onsider the du s t extin ction. The dust 
model is similar to lYaiima et al.l |201l3)) 6.g., rrid = 



0.01 X mg 



with a silicate- like mass density of 3 gem 



where m^, rrig and Z are dust mass, gas mass and 
metallicity. The size distribution is nd(ad) oc a^'^'^ 
IjMathis . Rumpl fc Nordsieclj Il977l ) with the dust radius 
of 0.1 — 1.0 \im. Since the assumed range of dust size is 
larger than the wavelength of Lyman limit, we assume the 
absorption efficiency Qiv) — 1 for the ionizing photons 
l|Draine fc Leelll984f) . 



2.2.2 UV Background Radiation 

The UVB can ionize the baryonic gas in galaxies, and then 
gas can be heated up to ~ lO** K. It would be ideal to com- 
pute the RT of both UVB and stellar radiation as the simu- 
lation runs, however in practice it is a very expensive calcu- 
lation. In general, the calculation amount of UVB transfer is 
large r than that of stellar radiation by 2 — 4 orders of magni- 
tude (jYaiima et al.|[2011al l . Therefore we fir st use the simple, 
four d ifferent models of UVB, similarly to iNagamine et al.l 
(|2010l ). We then compute the actual UVB transfer for some 
galaxies to check the validity of the UVB models. 

Our fi ducial simulations inclu de a uniform UVB with 
a modified iHaardt fc Madaul (|l996t ) (hereafter HM96) spec- 
trum (see 
place at z 



Dave et al.lll999l ). where the reionization takes 



6 as s uggested by the quasar observations (e.g.. 



iBecker et a"l]|200ll') and stellar radiative transfer calculations 
fe.g.. lSokasian et al.|[200^ ). 

We use the following four simulations with different 
treatment of UVB (but with the same initial condition) to 
examine the effects of UVB: 

(i) Optically-Thin Model: A uniform UVB radiation 
field with an optically thin approximation. 

(u) HMO. 5 Model (HM96 x 0.5): The ISM is optically 
thin to the same UVB, however the intensity of UVB is 
reduced to the half of the Optically- Thin model. 

(ifi) OTUV Model (Optically Thick UV): The ISM 
is optically thin to the UVB in the lower density regions 
with riH < O.Olpth = 6.34 x lO"'* cm~^, but completely 



optically thick in higher density regions of nn ^ O.Olpth, 
where pth is the threshold density, above which star 
formation is allow e d. Th e value of pth was determined by 
I Choi fc Nagamind (|2010l ) based on the observed SF cut-off 
column density of the Kennicutt law. The OTUV model 
implicitly assumes that the UVB cannot penetrate into 
the high density regions by self-shieldi ng. This type o f 
self-s hielding model was first propose d by Katz et aU (|l996l ) 
and iHaehnelt. Steinmetz fc RauchI (Il998ll to examine 

sr|r 



f(N-m) and DLA kinematics. Later INagamine et al.l (|2010l ) 
successfully reproduced the observed /(A^hi) by incorpo- 
rating t he OTUV model in t o cosmological si mulat ions. 
See also IPontzen et al] (|2008h : iFumagalli et al.l (|201lh and 
lAltav et all (|201ll '). 



(iv) No-UVB Model: UVB does not exist at all. 

Since we use the same initial condition for these runs, we 
can compare the effect of radiation halo by halo basis. After 
discussing the above four models in the beginning of Sec- 
tion [37l] we then move on to discuss the effects of radiative 
transfer in Sections I3.1.1l and l3.2.2l We compute the transfer 
of local stellar radiation in each halo as post-processing, and 
estimate the ionization structure in the equilibrium state by 
evaluating 



TiHi + r ' rlHi 



(1) 



where F^yg and T^^^^ are the photoionization rate by UVB 
and stellar radiation respectively; F*^ is the coUisional ion- 
ization rate; nc,nm and nun are the number density of free- 
electron, neutral and ionized hydrogen, respectively; ub is 
the total recombination coefficient to all bound excitation 
levels. The value of Fj.j^,. is estimated by the full RT calcu- 
lation, but Fyyg is computed with the above UVB model. 



3 RESULTS 

We first compare our simulation results to the re- 
cent observations of Hi c olu mn density functi on (e.g. . 



Prochaska fc Wolfe 
Prochaska fc Wolfe 



I2OO9I : iNoterdaeme et al.l l2009l : 

and examine the effect of 
UVB and stellar radiation on f{Nm). We then examine 
the Hi cross section of DLAs and LLSs, and study how 
they change with the UVB models and the RT calculation. 
Finally, we study the typical physical properties of DLA 
and LLS gas, such as SFR and metallicity. 

3.1 Column density distribution function 

The column density distribution function f{Nui,X{z)) is 
defined such that f{Nm, X)dNmdX is the number of ab- 
sorbers per sight line with H I column density in the interval 
[A^Hi, A^Hi + dA^Hi], and absorption distance in the interval 
[X,X -\- dX]. The absorption distance X{z) is given by 



X{z) 



/I I '\2 Ho 



dz'. 



(2) 



In practice, if the comoving box-size of the simulation is AL, 
then the corresponding absorption distance per sight-line is 



AX : 



-Ho 
c 



{l + zfAL. 



(3) 
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Figure 1. Column density distribution function at 2 = 3. The 
purple dash-long dashed line indicates the /{Nm) considering 
only star-forming galaxies that have young stellar sources (fage 
10^ yr). The orange dot-dash line indicates the /(Nm) consider- 
ing only massive galaxie s of ^ 5 X 10^'^ Mq. The obscrva- 



tion data points are from lPeroux et ah 1 1I2OOII. blue filled circles,) 
IPeroux et all l|2005l . black filled squares) O'Meara et alj l|2007) . 
magenta open circl e s). iProchaska fc Wolfel i200g|. green triangles) , 
iNoterdaeme et alJ l|200g| . cyan 

1I2OIOI. yellow shade) . The figure is similar to Figure 3 in 
iNagamine et aL I1 I2OIOI) . but the A'^HI range is extended to a lower 
value of Nm = 10^^ cm"'^ to discuss LLS. 



For example, for AL — 10 h ^ Mpc and 2: = 3, we have 
AX = 0.0534. 

Since the simulation box size is very small, we can safely 
assume that the DLA clouds do not overlap along a sight- 
line, and we can compute the A^hi distribution function by 
counting the number of projected grid-cells with column 
densities in the range [Nui, Nm + dNm] for all halos in the 
simulation. We use the same grid that we set around each 
dark matter halo for the RT calculation, with the grid cell 
size equal to the comoving gravitational softening length. 

We calculate the H I density in each grid-cell of volume 
(e/[l -f 2])^ after the RT calculation, and estimate Nui by 
projecting the H I mass distribution along the z-axis direc- 
tion as 



mp(l -I- z) 



(4) 



where e is the comoving gravitational softening length, rrip is 
the proton mass, and z is the redshift. The resulting /(Nni) 
is insensitive to the direction of projection. By this method, 
we are treating each projected grid-cell element as one line- 
of-sight. We have also checked that we are not missing any 
LLS columns in the outskirts of the halo by doubling the 
size of the grid for some of the halos. 



Figure[T]shows ,f{Nm) in the range of 17 < log TVhi < 22 
at 2 = 3 together with some observational data. The no- 
UVB run overpredicts the observed data points at all Nm 
range. In particular, the difference becomes as large as an 
order of magnitude at Nm ~ lO^^cm"'^. This shows that 
the ionization effect by the UVB is crucial in reproduc- 
ing the observed LLS and DLA column density distribution 
functions. On the other hand, the f{Nm) of the Optically- 
Thin and the HMO. 5 runs are close to the observed data 
at Ahi < 10^^'^ cm~^, however, falls below the observation 
data at 10^^'^ cm"^ < Ahi < 10^^ cm~^ by a significant 
amount. 

The OTUV run agrees with the observed /(Ahi) much 
better than the other models in a wide range of 10^^ cm^^ < 
Ahi ^ lO^'^ cm~'^. At lower column densities of Ahi ^ 
10^* cm"^, the OTUV run is lower than the OpticaU-Thin 
run, but it still fits with the observed range indicated by the 
yellow box. Thus, as far as the column density range of LLS 
and DLA is concerned, we conclude that the UVB treatment 
of the OTUV run successfully reproduces the observational 
data in a wide range of lO" cm~^ < Ahi < 10^^ cm"^. We 
note that the result o f the OTUV run in Fig ure [T] is slightly 
different from that in INagamine et al] (|2010l ). but this is be- 
cause the galactic wind model was slightly different between 
the two works. We hav e also checked the resu l ts wit h the 
updated UVB model of iFaucher-Giguere et al] (|2009l ). and 
found that the results are not very different. 

One notable feature in Figure[T]is that the OTUV run is 
higher than the observed /(Ahi) at Ahi > 10^^ ''. These high 
column density sight-lines mainly come from haloes of higher 
dust-to-gas ratio. In such halos, the molecular hydrogen H2 
is cre ated efficiently, and t he H2 mass fraction can be > 0.1 
(e.g.. lHirashita et al.ll2005l ). Because our current simulations 
do not include H2 formation, we may overpredict /(Ahi) at 
high column densities of Ahi ^ 10^^'^ cm~^. 

We find that the stellar radiation transfer does not af- 
fect the shape of /(Ahi) (overlapping blue dashed line and 
red solid line in Figure [1]), whereas the UVB model changes 
the /(Ahi) significantly. We will discuss the effect of stel- 
lar radiation on DLAs and LLSs in more detail later in the 
paper. 

The change of column density by UVB in the OTUV 
ru n is similar to those fou nd by previous works (e.g., Fig. Al 
m iFumagalli et al.ll201lh . On the other hand, the effect by 
st ellar radiation in our sim ulation is much smaller than that 
of IFumagalli et all (|20lj ). There could be several possible 
reasons for this difference as follows. In our simulations, the 
high-density gas clumps near young stars can effectively ab- 
sorb stellar radiation, and block the propaga tion of stellar 
radiation over extended regions in a galaxy (|Yaiima et al.l 
l2011al ). Another possible reason could be different numer- 
ical methods: our RT method is based on a ray-tracing 
method, whereas Fumagalli's method is Monte Carlo. The 
ray-tracing method can estimate the ionisation structure 
by point sources more accurately, while the Monte Carlo 
method can cause s ome artificial inh omogeneity (see the 
comparison work by llliev et aL I I2OO6I ). Finally, the spatial 
resolution is always a concern, and we note that the resolu- 
tion of Fumagalli's AMR simulation in high-density regions 
is better than ours. Unfortunately it is very difficult to tell 
the exact resolution effects unless we compare the two codes 
closely side by side. 
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Figure 3. Panel (a) : Neutral fraction of hydrogen gas as a function of number density for the same halos as in Figure [2] The black 
dashed lines and blue solid lines indicate the equilibrium solution of ionization balance between the UVB ionization and recombination 
in the Optically-Thin model and the OTUV model, respectively. The red filled circles are obtained by including direct radiative transfer 
calculations of UVB, stellar radiation and coUisional ionization. The blue filled circles of Halo C are the neutral fraction with UVB and 
coUisional ionization only (i.e. without stellar radiation transfer). The green filled circles are the neutral fraction with only UVB (i.e. 
without stellar radiation and coUisional ionization). Panel (b) : Neutral fraction of gas clumps irradiated by the radiation from young 
star clusters at the indicated distance £ as a function of hydrogen number density. See text for the details of the calculation. Different 
colored lines indicate different conditions as indicated in the legend, where N-py^ [s~^] is the emissivity of ionizing photons from star 
clusters, and I [kpc] is the distance between the star clusters and the gas clump. 



Another major source of different results might be the 
different mass range of simulated galaxies in each simulation. 
iFumagalli et al.l (|201ll ) focused on only seven star-forming 
galaxies with Mh > 5 x Mq . On the other hand, our sim- 
ulations cover a wider range of galaxy masses, and there are 
numerous low-mass, faint galaxies without young star clus- 
ters that contribute to the DLA cross section. In those low- 
mass haloes, the local stellar radiation does not impact the 
ionization of gas, and Fumagalli et al.'s simulation did not 
include these low-mass halos. In our simulations, about half 
of the total DLA cross section is contributed by faint low- 
mass haloes where there is no young stars in our simulations. 
The DLAs in these faint haloes are unaffected by the stellar 
radiation. This could explain the little effect of local stellar 
radiation on f{Nm) in our simulations, and more effect in 
Fumagalli's. In addition, Fumagalli et al.'s f{Nm) showed 
good agreement with observations at A'^hi ^ 10^"'^ cm"'^, 
but underpredicted at 10^** < A'hi < 10^° cm"^. They sug- 
gested that the reason for the underestimate might be their 
limited galaxy sample with Mh ^ 5 x 10^° Mq. For compar- 
ison, here we investigate the effect of galaxy mass limit in 
our simulation in Figure [1] When we consider only the mas- 
sive galaxies in haloes with Mh > 5 x lO^^M©, then /(Ahi) 
decreases for DLAs with a similar slope to the all-halo case, 
and f{Nm) is significantly lower and more flattened for the 
LLS. Integrating f{Nm) over the column density range of 
LLS, we find that the contribution from the massive galax- 
ies to the line density of LLSs is only '^20% of the total. 
Hence, we confirm that the LLSs mainly come from lower- 



mass galaxies, and t he discrepancy o f ,f(Nm ) with observa- 
tion at LLS range in IFumagalli et ahl (|201lh is probably due 
to the mass limit of their galaxy sample. 

In addition, we compute the f{Nm) using only the star- 
forming galaxies in our simulations, and find that the result 
is in-between the all-halo case and only massive galaxy case. 
This result suggest that LLSs mainly originate from lower 
mass galaxies in haloes with Mh < 5 x 10^°Mq and lower 
star formation rate. We will discuss this point more quanti- 
tatively in Section [3.61 

3.1.1 Radiation field by UVB transfer 

Next we perform direct radiation transfer calculation of 
UVB to validate the results of the OTUV run. We perform 
this calculations by irradiating UVB from the boundary sur- 
faces of calculation boxes with rig (~ 10®^^") rays, and also 
consider the stellar radiation transfer simultaneously. 

Figure [2] shows the photoionization rate versus gas den- 
sity in each grid cell from direct UVB transfer calculation. 
Here we pick four halos randomly from the halo sample of 
Mh < 10^^ Mq in the OTUV run, and recompute the trans- 
fer of UVB directly in each cell as shown by the red points. 
The higher mass haloes are excluded from our sample, be- 
cause the UVB RT calculation takes too long to process 
such a large grid (>100^) with current computers. The red 
points in the figure are derived from the RT calculations 
including both the UVB and stellar radiation. The black 
dashed and blue solid lines are the photoionization rates 
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Figure 4. Physical cross section of DLAs {panel [a]) and LLSs {panel [b]) as a function of halo mass at z = 3. Each point represents 
a halo with star-forming galaxies, including the stellar RT. Filled triangles indicate the median values in each mass bin. In the bottom 
right panel, the results from different runs are compared: Optically- Thin (red), OTUV (blue), HMO. 5 (green) and No-UVB (magenta) 
runs, respectively. 
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Figure 2. Photoionization rate by UVB radiation (Fuvb) as a 
function of hydrogen number density for four different halos in 
the OTUV run (Halo A : = 9.7 X 10^°Mq, B : = 2.5 x 
101°Mq, C : Mh = 1.5 X 101°Mq, and D : = 9.4 X 10^ Mq). 
The halos are taken from the OTUV run, but we are solving the 
transfer of UVB radiation directly in each cell as shown by the red 
points. The black dashed lines and the blue solid lines indicate 
the ionization rate in the Optically-Thin model and the OTUV 
model, respectively. 



by the Optically-Thin model and the OTUV model, respec- 
tively, obtained from the equilibrium calculation. Note that 
the photoionization rate in this figure is of UVB radiation, 
r^yg in equation (1), not including stellar radiation FJ-j^j.. 

We find that the Fuvb has a relatively sharp tran- 
sition from a high plateau at low densities to a much 
lower value at logriH ~ —2.2, and the self-shielding thresh- 
old density (nj{j^ = 6 x 10~^cm~^, indicated by the blue 
solid line) adopted in the OTUV model correctly cap- 
tures this sharp transition . The value of nj{j^ was deter- 
mined by iNagamine et alj (|2010l ) by examining the agree- 
ment with the observed data on f{Nm), and one of the 
main purposes of this paper is to validate this choice of 
n^y^ through direct RT calculations as demonstrated here. 
INagamine et al.l (|2010[ ) also gave some justification for this 
particular choice of n^^, referring to prior radiative trans- 
fer simulations. In addition, the adopted threshold den sity is 
close to that of iHaehnelt. Steinmetz fc RauchI (|l998h . who 
suggested the density threshold model based on the tight re- 
lation between column density and density. The region with 
riH ^ lO"'^ cm~^ mainly contributes to the column density 
of A^Hi > 10^'' cm~^, which can self-shield from UVB. 

Some limited number of cells at tihi > 10~^ cm"'' for 
Halos A & C show relatively high values of F^, deviating 
from the main photoionized branch of gas. These are cells 
that are probably ionized by the local stellar radiation, and 
therefore allows UVB to penetrate into deeper densities. In 
some regions near star-clusters, the stellar radiation can in- 
crease the photo-ionization rate than UVB alone by orders 
of magnitudes, and ionize star-forming regions. However, the 
impact is not so significant for a statistical quantity such as 
f{Nm) (Figure[l|. 

Figure [3^ shows the neutral fraction of hydrogen gas 
(Xhi) in each cell as a function of number density for the 
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same four halos shown in Figure (2] In Halos A and C, 
the stellar radiation boosts the ionization degree than the 
optically-thin equilibrium values, populating the red points 
below the black dashed curve. In Halo B and D, the young 
star clusters are distributed in high-density peaks, hence 
most ionizing photons are absorbed on the spot. In addition, 
the neutral fraction Xhi in low density region is approxi- 
mately fitted by Xhi ~ aeWH/r^yg, therefore the neutral 
fraction is simply proportional to the number density of hy- 
drogen. In addtion, we calculate the neutral fraction with 
different ionization processes (green : UVB only, magenta : 
UVB -I- coUisional ionization, and red: UVB + stellar ra- 
diation -I- collisional ionization). It shows that the collision 
boosts the ionization degree at ~ 10^* cm^'^ (magenta 
points extending below the dashed curve) , and the stellar ra- 
diation boosts the ionization degree even at higher density 
regions, as signified by the wide distribution of red points 
below the dashed curve. 

We can further understand the transition from ionized 
to neutral state by solving for the balance between pho- 
toionization by stars and recombination in the one-zone ap- 
proximation, picturing the gas clumps near young star clus- 
ters. In Figure Is]], we estimate the photoionization rate by 
Tatar = A'phcroexp(— r)/47rZ^ , where TVph is the emissivity of 
ionizing photons from star clusters, ctq is the ionization cross 
section at the Lyman limit, / is the distance between star 
clusters and gas clumps. The optical depth r is roughly es- 
timated by r = aonml, where nui is the number density of 
neutral hydrogen of the gas clump. The A^ph of a star cluster 
in our simulations can be ~ 10^""^^ s~^, and the distance 
I can be proper ~ 0.5 — 1.0 kpc, corresponding to the cell 
size of RT calculation in our simulation. In Figure [SJj, the 
neutral fraction steeply rises at the density where r becomes 
~ 1. With increasing photon emissivity, the threshold den- 
sity moves to a higher density as the radiation is able to 
sink into higher densities, and the ionization degree at lower 
densities become stronger. With the above parameter val- 
ues, we find that the star clusters cannot ionize gas clumps 
with TiHi S> 1 cm~'^. In our simulations, the gas density near 
star clusters can be > 1 cm~"^ in massive halos and block the 
ionizing radiation from stars. From these analyses of UVB 
radiation transfer, we confirm that the stellar radiation does 
not strongly affect the neutral fraction of high-density gas 
at riH > TT^th^, and the UVB treatment of the OTUV run is 
reasonable. 

3.2 Cross Sections of DLAs and LLSs 

Once the Hi column density of each cell in the projected 
plane is obtained, we can estimate the DLA cross-section 
of each halo by counting the number of grid-cells that ex- 
ceed A'hi = 2 X lO^'' cm~^ and multiplying by the unit area 
(e/(l + z))'^ of the grid-cell. In this paper, we discuss the 
cross section in physical units. 

First, we examine the effect of different UVB models on 
the cross section. Figure [4^,b show ctdla and (tlls of each 
halo as a function of halo mass for different UVB runs. Each 
point in the figure represent a halo with star-forming galax- 
ies, including the effect of stellar RT. Each run has ~ 2000 
galaxies in the simulation box at 2: = 3, and about 200 of 
them have young star clusters (tagc s; 10^ yr) that can af- 
fect the hydrogen ionization (hereafter we call these galax;ies 



Run 


a(DLA) 


a (LLS) 


13 (DLA) 


13 (LLS) 


Optically- Thin 


0.65 


0.61 


2.12 


3.59 


HMO. 5 


0.67 


0.56 


2.23 


3.66 


OTUV 


0.70 


0.53 


3.06 


3.41 


OTUV (no-star) 


0.72 


0.55 


3.11 


3.46 


no-UVB 


0.66 


0.30 


3.48 


3.94 



Table 1. Best-fit parameters of tlie fitting function (Eq. [5]) for 
'''DLA (Tlls fo'' different UVB models. The value of a is the 
slope of the power law, and f3 is the normalization point at Mjj = 
10^2 Mp). 



"star-forming galaxies"). Here, we select the galaxies with 
young stars from each run, and compute the RT of local 
stellar sources. 

We find that the cross section increases with increasing 
halo mass for both DLAs and LLSs in all UVB models. For 
a quantitative discussion, we plot the median value in each 
bin with a triangle symbol, and then fit a power law ctdla oc 
M^, assuming a functional form of 

logcr = a(logMh - 12) (5) 

where Mh is the dark matter halo mass in units of solar 
masses, and the values of slope 'a' and normalization are 
determined by least-square fitting. We use the cross section 
at Mh = lO^^M© as the anchor point of /3. The best-fit 
values of a and /3 of all runs are summarized in Table [T] 

The bottom right panels of Figure |4^,b compare the re- 
sults of different UVB models. For both DLAs and LLSs, the 
cross section is the greatest for the No-UVB run, and it drops 
dramatically once the UVB is turned on. The difference be- 
tween the Optically-Thin run and the HMO. 5 run is small, 
suggesting that the gas is ionized to relatively high densities 
under the optically-thin approximation, irrespective of the 
specific strength of the UVB. For DLAs, the OTUV run has 
much higher cross sections than in the HMO. 5 run, because 
the high-density clouds with n > are the main contrib- 
utor to ctdla, and the OTUV model increases the neutral 
fraction of such gas over the optically-thin approximation. 
In fact, the value of P for the OTUV run is much higher than 
that of HMO. 5 run, but the slope is similar in all runs with 
a « 0.65 — 0.72 for DLAs. This suggests that ctdla is af- 
fected proportionally at all halo mass range by the different 
UVB models. Obviously t he result of the Optically -Thin run 
is very similar to that of iNagamine et al.l (|2004bl '). as they 
both used the optically-thin approximation in similar simu- 
lations (the main difference is the galactic wind model). 

The biggest difference between (tdla and (Jlls is that 
the OTUV run (blue line) has almost the same ctlls as the 
Optically- Thin run (red line) , whereas (Jdla are very differ- 
ent in the two runs. This implies that the high-density gas 
with n > nYi^ is not responsible for the LLS. The HMO. 5 
run has a slightly higher (Jlls than the Optically- Thin run 
as expected. We also find that the slope a is systematically 
shallower for the LLSs compared to that of DLAs, and the 
normalization /3 is higher for the LLSs. This is expected, 
because LLSs are more extended in the outskirts of halos. 

If we look at the distribution of data points in Fig- 
ure|4^,b more carefully, one can see that there is a slight cur- 
vature to the median points, and a straight power-law does 
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Figure 5. Comparison of DLA cross section with other theoret- 
ical works. The black dotted and dashed lines are t he resu lts of 
P3 and Q3 models, respectively, bv lNagamine et al. I ll2004al) . The 
black long-dashed and dot-dashed lines are the strong wind (SW) 
and weak win d (WW) runs , resp ectively, in the GADGET SPH 
simulations of lTescari et al.l l|200gf l . The black filled circles arc the 
cr oss sections of seven si mulated galaxies in the AMR simulation 
of lFumagalli et~al The cya n and orange shade en close the 

data points from the simulations of lPontzen et all | |2008| . Gasoline 
SPH simulation) and lCeiJ l|201Cll . Enzo AMR simulation), respec- 
tively. 



not necessarily provide a perfect fit. At low halo masses, the 
median points fall below the power-law fit. At the interme- 
diate to high masses, the median points are slightly higher 
than the power-law fit, and the cross section seems to flat- 
ten out at the highest masses at Mh 



10^^ A similar 



flattening has been observe d in other works (jPontzen et al.l 
I2OO8I : iFumagaUi erallboill ). as we will discuss in the next 
section in more detail. 



3.2.1 Companson with Other Simulation Results 

In Figure (5] we compare the results of DLA cross 
section from several authors using different simulations 
llNagamine et al. l2004bl: IPontzen et aljlioosl : iTescari et al] 
I2OO9I : ICenlbOld : [FumagaUi et aLlbOHl ). Although there are 
small differences in the slope and normalization, most of the 
results are bracketed by our no-UVB run and Optically- Thin 
run. There a re significant o v erlap b etween the result s of ou r 



OTUV run, IPontzen et all (j2008l ). iFumagalli eTal 



2011), 



20091), 



and the WW (weak wind) model of iTescari et al 
whi ch is very encouraging . 

iPont zen et ahl (|2008l ) considered a crude self-shielding 
model, so it is natural that their re sult is close to the OTUV 
run. Although ITescari et al.l (|2009l ) also set the self-shielding 
threshold density criteria above which UVB cannot pene- 
trate, their nj^^ is about one magnitude higher than our 
value. Hence, their SW (strong wind) model should give 



I I I I I I I I I I I I I I I I I 

m 




LLS 



9 10 11 12 
Log M, [MJ 



1 

0.8 

°^ 0.6 

0.4 

0.2 


1 

„0.8 

o in 

^ 0.6 

0.2 


9 10 11 12 
Log M, [MJ 

Figure 6. Upper panels : The ratio of cross sections of DLAs 
and LLSs in each halo with (cr) and without {u^') stellar RT: 
i.e., ct/u^. Lower panel : The mass ratio of neutral hydrogen in 
each halo with (Mjji) and without (M^j) stellar radiation, i.e., 
Mhi/M°j. The red points are the mean value in each bin with 
the bin size of 0.5 dcx. 



lower o-DLA than our OTUV run. The SW and WW model 
have constant wind particle velocities of 600 and lOOkms"^ 
as a stellar feedback, respe ctively. The d i fferenc e between 



the SW and WW mod el oflTescari et ahl (|2009|) 



tent with the results of lNagamine et alj (|20(4bl) 



Using an Enzo AMR simulation, ICenI (j201Gl ) reported 
a much higher (tdla than other simulations for halos with 
Mh = 10^° - 10" Mq, and the reason for this difference is 
not very clear. Some of the possibilities are: 1) their star 
formation model is considerably inefficient in converting the 
neutral gas into stars, 2) their self-shielding treatment is 
allowing too much neutral gas in massive halos, and/or 3) 
their stellar feedback is inefficient and does not heat up the 
ambient neutral gas. They used only two zoom-in simula- 
tions (a high-density cluster region and a low-density void 
region) to bracket the mean density result, therefore they 
only bracketed the observed f{Nm) without reproducing 
the normalization of f{Nm). It appears that his result is 
more consistent with that of our No-UV run rather than the 
OTUV run. 

We find that the Fumagalli's fit to the ctdla for low- 
mass halos is much steeper than other simulations, which 
could be due to a small galaxy sample in their zoom-in sim- 
ulations (only seven galaxies). We discussed the effect of 
limited galaxy sample already in 13.11 and Figure [T] 



3.2.2 Effect of Stellar RT on Hi Cross Sections 

The upper panels of Figure |6] compares the ratio of DLA 
and LLS cross section with and without stellar radiation. 
The data points are only for those halos that contain star- 
forming galaxies in the OTUV run at z = 3. We see that 
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Figure 7. Two dimensional map of A^hi for four halos in the OTUV run at 2 = 3. Each column shows different treatment of radiation: 
column (a): no-ionization; (fe): UVB + coUisional ionization, (c): UVB + star + coUisional ionization. The column (d) is the same as column 
(c), but for the same halo in the Optically-Thin run. Each row corresponds to one halo with a panel size ^phys in physical units: Halo 
(1): Mh = 6.7 X 10" Mg.^phys = 300 kpc; Halo (2) = 1 X lO" Mg.-^phys = 159 kpc; Halo (3) = 2.6 X 10^° Mg.fphys = 80 kpc; 
and Halo (4) Mjj = 3.3 X 10^ ^Q^f-piiys = 47 kpc. The white vertical tick mark in column (d) indicates a scale of proper 30 kpc. 



both (Tdla and (Tlls of lower mass haloes are reduced by the 
stellar RT, but the higher mass haloes are not affected very 
much. The slope 01 and normalization /3 for the power-law 
fit are very similar in the two cases of with and without the 
stehar RT. Since young stars distribute near high-density 
gas clumps, most of ionizing photons from stars are blocked 



by the gas clumps, and the region behind high-density gas 
clumps are not irradiated by the stellar radiation. As a re- 
sult, the stellar radiation cannot propagate over long dis- 
tances in the halo, and the stellar radiation does not affect 
the Hi cross section very much in high-mass halos. In ad- 
dition, although relative mass of neutral hydrogen gas to 
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stellar mass somewhat decreases as increasing halo mass, 
the amount of dense gas clouds around star-forming regions 
increases with halo mas s by frequent minor merging pro- 
cesses (e.g., Figure 1 in lYaiima et al.ll2011al ). The photon 
absorption process and the distribution of stars and gas are 
determined by complicated nonlinear process, which can be 
traced by cosmological hydrodynamics and radiative trans- 
fer simulations. 

The lower panel of Figure [S] shows the ratio of H I mass 
in each halo with {Mm) and without (M^i) the stellar RT, 
i.e., Mui/Mm- It shows that the lower mass halos lose more 
Hi masses than the higher mass halos, consistently with 
Figure |6l The H I mass is reduced by about 30 per cent by 
the stellar RT for halos with Mh < 10^" M©, while it is 
only about 7 per cent reduction for higher mass halos with 
Mh > 10" Mq. 

Some analytical works discussed the relative strength 
of local stellar radiatio n to UVB under the assumption o f 
simple gas distribution (|Mhalda-Escudell2005l : [Schavdl2006l ). 
iMiralda-Escudi (|2005l ) showed that the effect of stellar ra- 
diation is negligible compared to UVB for LLSs, and that 
it can be comparable for DLAs a s an u pper limit. On the 
other hand, the model bv lSchavd (j2006l ) indicated that the 
local stellar radiation can be dominant even for LLSs, as 
well as for DLAs. However, for inhomogeneous gas distri- 
bution, the effect of stellar radiation strongly depends on 
the structure of gas around stars, and it requires RT sim- 
ulations for proper assessment. The above analytical works 
did not discuss the direct impact of local stellar radiation 
on DLA/LLS cross sections and /(A^hi), and focused only 
on the relative strength between local stellar radiation and 
UVB. In our simulations, the ionizing flux from local stars 
can be much higher than the UVB similarly to the previ- 
ous analytic estimates. However, the neutral clouds survive 
from io nization by self-shield ing and shadowing effect (See 
Fig.l of lYaiima et~al]|2011al '). As a result, the local stellar 
radiation cannot change the ionization structure in galax- 
ies very much (particularly in massive galaxies), and hence 
the H I cross section does not change largely. On the other 
hand, the UVB enter from outside the halo into low-density 
region, therefore they can irradiate over greater volume. 

3.3 Structure of DLAs and LLSs 

As we see in Figure [U the UVB can drastically change the 
rate-of-incidence of DLAs and LLSs. In this subsection, we 
investigate how the UVB and stellar radiation change the 
spatial distribution of A^hi. Figure [7] shows the two dimen- 
sional map of Ahi distribution for four halos in the OTUV 
run at 2 = 3. The column (a) is without any ionization by 
UVB and stellar radiation, corresponding to the No-UV run. 
The column (6) is with the UVB RT and coUisional ioniza- 
tion, but no stellar RT. The column (c) is with the stellar 
RT in addition. The column (d) is the same as column (c), 
but for the same halo in the Optically- Thin run. 

Comparing columns (a) and (6), we see that the UVB 
changes the Ahi distribution drastically. Much of the LLSs 
in the outer region are ionized by the UVB and becomes 
lower than Ahi ~ lO^^cm"^. The higher columns with 
Ahi lO^'''^ survive from photo-ionization by UVB because 
of the high recombination rate in high density regions, but 
the DLA cross section is significantly reduced. 
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Figure 8. Cumulative probability distribution function (PDF) 
of DLAs and LLSs as a function of separation distance (i.e. impact 
parameter) between DLAs (LLSs) and center-of-mass of the near- 
est galaxy in projected 2D plane in physical unit. The colored lines 
are our simulation results in this work. The black dash and dot- 
dash lines arc PDFs with the constant- velocity ( CW) wind model 
-I- op tically-thin UVB approxima tion (P3 run inlNaeamine et al.l 
\20Qlh and the CW-I-OTUV run llNagamine et al.ll2O10l) 7 respec- 
tively. 



The DLAs correspond to the yellow and red regions in 
Figure [T] The DLA distribution does not change very much 
from column (&) to (c), and it is not affected by the stellar 
RT very much as we discussed in § 13.2.21 

Comparison of columns (c) and (d) reveals the effect of 
optically-thin approximation of UVB in the Optically- Thin 
run relative to the self-shielding effect in the OTUV run. 
It shows that the UVB photo-ionizes the gas too deeply 
in the optically-thin approximation, and the self-shielding 
model increases the DLA cross section significantly. From 
the discussion in § 13.11 we argue that column (c) is more 
realistic than column (d) based on the comparison of /(Ahi). 

In the OTUV run (column [c]), DLAs are extended over 
~ 10 kpc scales in all halos, which is a stark contrast to the 
very compact appearance in the Optically- Thin run (column 
[d] ) , a nd the LLSs are surround ing DLAs on ~ 30 — 60 kpc 
scales. iNagamine et all (|2004bh showed the DLA distribu- 
tion similar to the column (d) in their SPH simulations with 
optically-thin approximation and argued that DLA distribu- 
tion is very clumpy. However based on the present work, it 
seems that DLAs are more broadly and smoothly distributed 
around LBGs at 2 = 3 on ~ 10 kpc scales. The highest col- 
umn density systems with A^hi > 10^^ cm~^ are preserved 
well even in the Optically-Thin run, as we saw in Figure [T] 

To quantify the relative distribution of DLAs, LLSs and 
galaxies in different UVB models, we compute the impact 
parameter distribution of DLAs/LLSs relative to galaxies in 
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each run. Figure [8] shows the cumulative probability distri- 
bution function (PDF) of DLAs and LLSs as a function of 
separation distance (i.e. impact parameter) between DLAs 
(LLSs) and the center-of-mass of nearest galaxy in projected 
2D plane in physical unit. The colored lines are from the sim- 
ulations used in this paper. The OTUV run shows a much 
shallower PDF than the Optically-Thin run. About fifty per 
cent of DLAs in the OTUV run distribute within ^ 7 kpc 
from galaxy center, while those in the Optically- Thin run 
are concentrated within ~ 3 kpc. This is because the UVB 
cannot penetrate deeper into the galaxy center due to self- 
shielding, and hydrogen gas near galaxy center can survive 
as DLAs. Of course, no- UVB run has the widest PDF due to 
the overproduction of DLAs as noted earlier in this paper. 

In addition, we compare the results with different galac- 
tic wind models. As described in Section 2.1, we use the 
" Multi-component Va r iable Velocity" (MVV) wind model 
of IChoi fc Nagamind (|201ll ) in the current simulations. 
Some previous works used the constant-velocity wind (CW) 
model. The black dashed and solid lines are the PDFs 
from the CW model -I- optically-t hin UVB approximation 



fPS run in iNagamine et all |2007| ) and CW-hOTUV run 
(jNagamine et al.l 20101 ). respectivelv. We find that the DLAs 



in the CW model are concentrated near the galaxy center 
more than the MVV model. In the CW model, the lower 
mass galaxies have higher win d velocities than in the MVV 
run jChoi fc Nagamind l201lf ). therefore the DLAs abun- 
dance is somewhat suppressed. In such a situation, only 
high-density gas around galaxy center can survive as DLAs, 
resulting in a steeper PDF. In the case of LLSs, the PDFs 
are in general much wider than those of DLAs'. 



3.4 Stellar Mass and Metallicity Distribution 

We show the projected maps of stellar mass, A^hi, and metal- 
licity in Figure [51 We find a strong correlation between high- 
A^Hi systems and stellar mass distribution as expected. The 
metallicity map is obtained by summing up all metal mass 
along each line-of-sight and dividing by the total gas mass 
in the same sight-line. The high density region is effectively 
enriched by the metals created by SNe. The massive halos 
(1) and (2) show high metallicities around the high density 
regions and DLAs. 

Galactic wind can also carry high-metallicity gas out- 
side the galaxies. All bottom panels of Figure [9] show that 
the near-solar metallicity gas is carried well outside of the 
stellar distribution, up to ~ 30 — 50 kpc away from stars. The 
effect of wind appears stronger in halos (3) and (4), but this 
is simply due to the smaller physical scale of the panels for 
halos (3) & (4) compared to (1) & (2). Careful examination 
shows that halos (1) & (2) also show a similar degree of cir- 
cumgalactic enrichment as halos (3) & (4). For example, for 
halo (1) there is a bubble of sub-Z0 gas extending toward 
north-east direction, about 60 kpc away from the center of 
galaxy on the right-hand-side. For halo (2), a similar fea- 
ture can be seen for a near-Z© gas extending out towards 
north-east direction, again about 60 kpc away from the cen- 
ter of galaxy on the right-hand-side. For halo (3), there is a 
large bubble of near-Z© gas extending towards south-west 
direction, about 30 kpc away from the stars. For halo (4), 
a similar feature can be seen towards north direction. For 



halos (1) — (3), notice that the direction of wind propagation 
is nearly perpendicular to the chain of galaxy distribution. 

Our MVV wind model reproduces the observed scaling 
between wind velocity and galaxy SFR: i;^, oc SFR^^^ oc 

1/3 

Mi, , which allows greater wind velocity for more massive 
galaxies. This is compensated by the deeper potential well of 
more massive galaxies, therefore the enriched metals reach 
similar distances for all halos as we discussed above. 

The bottom two rows of Figure [9] also show the virial 
radius of each halo in white circles. One can see that the 
DLAs are about the same or slightly smaller than the virial 
radius, whereas the metal enrichment clearly goes beyond 
the virial radius. For halos (1) & (2), there are two major 
galaxies in these massive halos, and each of the galaxies is 
a host of DLAs that extends to the scale of virial radius. 



3.5 Cumulative DLA Rate of Incidence 

Using the fitting function of DLA cross section derived in 
§ 13.21 we are now able to estimate the cumulative num- 
ber of DLAs per unit redshift. However, there is the mass- 
resolution limit in a simulation, and we cannot resolve the 
lower mass haloes that may be DLAs host. It may cause 
the underestimation of the DLA abundance. In addition, it 
is hard to produce very massive halos in a limited simula- 
tion volume. To overcome this resolution limit, the DLA 
abundance has been estimated by combining a theoreti- 
cal fit to the dark matter halo mass function and the de- 
rived relationship betw e en DLA cross-section and halo mass 
jGardner et al.lll997al lbl l200ll : iNagamine et al.ll2004d . l2007l : 



iTescari et al.ll2009h . Here, we follow the same method. The 
cumulative number of DLAs is estimated by 

dNuLA^^ M,2) = ^ f°° ni,{M',z)auLA{M',z)dM' 



dz 



(6) 

where ?ih(Af ,z) is the dark matter h alo mass funct i on, fo r 
which we use the function derived in [Tinker et al. 1 (|2008D . 
and dr/dz = c///o-\/r2m(l + z)'-^ + ^a- 

Figure [TOl shows the cumulative DLA number as a func- 
tion halo mass for different UVB models. The yellow shaded 
region is the observational estimate log(dn/dz) = — 0.6 ±0.1 
at z = 3 from SDSS data of lProchaska et al.l (|2005h . We find 
that the cumulative abundance largely depends on the UVB 
models. The OTUV run agrees well with the observed data 
, if the cross section steeply drops off at Mh < 10^ Mq . In 
this work, we set the lower limit of halo mass to lO'' Mq for 
the above integration, because the DLA cro ss section steeply 
declines at halo masses of 10*-^ - 10^ Mq (INagamine et all 
l2004al : |Pont zen et al ] |2008l ). However, there is some ambigu- 
ity in the DLA formation in lower mass haloes, and we would 
overpredict the abundance if we were to integrate down to 
10*'^ Mq for the above integration. 

On the other hand, the No-UVB run clearly over- 
produces DLAs. The Optically- Thin and HMO. 5 runs un- 
derpredict the observation by some factor. The result 
of the Optically-Thin r un is slightly lower than that of 
INagamine et al.l (|2004bl . [2OO7I '). The new MVV wind model 
used in this work shows the smaller ctdla a t higher halo 
mass. In addition, the halo number density by [Tinker et al.l 
(|2008h with WMAP 7-year parameters is slightly lower than 
the mass function used in their paper. Hence, it results in 
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Figure 9. Two dimensional map of projected stellar mass (top row), A^hi (middle row), and metallicity (bottom row) for the same four 
halos shown in Figure [7] in the OTUV run. Each column with the number index corresponds to the same halo in each row of Figure [T] 
and the panel sizes are the same as in Figure [71 The color scheme of stellar mass is normalized in each panel. The middle row is the 
same as the column (c) in Figure [T] but we show it again for comparison. The white bar in the middle row shows the proper 30 kpc. 
Note the different panel sizes for each column. The metallicity map is obtained by summing up all metal mass along each line-of-sight 
and dividing by the total gas mass in the same sight-line. The white circles indicate the virial radius of each halo: physical 68, 37, 23 and 
11 kp c for halos (l)-(4), respectively. The virial radius is estimated by Rvir — 144(1 + z^i^)-'^ {Mt^/W^'^ Mq)^/^ kpc (e.g.. lMo fc White! 
l2002t) . and here we assumed z^-„ = 



smaller DLA abundance even with the same optically-thin 
UVB prescription. Although our (Tdla in th e OTUV run 
is som ewhat larger than that of SW run in iTescari et al.l 
(|2009D . our result still agrees with the observation data well. 

3.6 Typical physical quantities of DLA host 
galaxies 

In this section, we examine the typical physical quantities of 
DLA host galaxies in the OTUV run. Figure [TT] shows the 
o"DLA-weighted probability distribution function (PDF) as a 
function of certain physical quantity. For example, PDF oc 
dcTDLA/rflog Afh for the top left panel. The histograms are 
normalized such that the area under the histogram is equal 
to unity. The dashed and dotted lines indicate the median 
and mean values weighted by the cross section, respectively. 
We find the median halo mass by rank ordering the halo 
mass according to the a;-axis quantity, and determining the 



halo at which the stacked cross section being half of the total. 
We see that the mean values are biased towards higher values 
and not representative of the entire distribution, except for 
the bottom right panel where the distribution is strongly 
peaked at a certain metallicity. Therefore median value is 
the better indicator of typical property of DLA hosts. 

The upper left panel of Figure [TT] shows that the me- 
dian halo mass of DLA hosts is Mh = 2.4 x IO^^Mq, and 
about half of the total (Tdla is contributed by the halos 
with Mh = 6.6 X 10^ - 1.1 x lO" Mq, centered around 
the median (±25% around median). Although (Jdla in- 
creases with increasing Mh, the number density of massive 
halos decreases at the same time. This effect pushes the 
median DLA halo ma ss to be in the medium mass range. 
iNagamine et al.l |20o3) also presented the median DLA halo 
mass for various wind models, and our results are similar to 
their P3 or Q3 model. The cyan shade indicates the mass 
range of DLAs host galaxies suggested by clustering analysis 
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Figure 10. Cumulative abundance of DLAs per unit redshift 
as a function of halo mass for the Optically-Thin (dotted line), 
HM0.5 (dot-dash hnc), OTUV (red solid hne) and No-UVB (dash 
line) runs at 2 = 3. The yello w shaded region shows t he observed 
cumulative DLA abundance of lProchaska et al.l ([200^ from SDSS 
data. 



ijCooke et al.ll2006t ). ICooke et al.l (|2006l ) showed that mass of 
galaxies hosting DLAs were 10^ '' < M < 10^^ Mq from 
DLA-LBG cross-correlation. Our simulation res ults agree 
well w ith Cooke's result, and the simulation work ()Lee et al.l 
In addition, we derive the bias parameter from the 
median mass. Here we use an ellipsoidal col lapse model 
ijSheth. Mo fc TormenllioOll : IMo fc Whitell2002l ) . which can 
reproduce results of A'"-body simulations well. We derive the 
rms fluctuation of mass density at a mass scale M, (j{M) 
by using the m edian mass of DLAs, and a ssign it in the 
equation (8) of ISheth. Mo fc TormenI (|200ll ). As a result, 
we obtain the bias of 6 = 1.69 (Table [2|. Thi s bias is also 
in the observed range of 1.3 < 6dla < 4 by I Cooke et al.l 
(|2006l ). 

The distribution is relatively flat as a function of SFR. 
The typical LBGs have SFR ^ I - 2Q Mq yr~\ and the 
fraction of ctdla contributed by such LBGs is ~ 30 per 
cent based on our simulation. The median SFR in the 
top right panel of Figure 1111 is 0.3Mmvr~^. Using the 
Equation (1) in iKennicuttI (| 19981 ) for the Salpeter IMF, 
this median SFR corresponds to ~ 0.02 L* of LBGs at 
2 = 3, which is much fainter than the observati onal flux 
limit (0.1 1/,) in the recent LBG survey of , Roddy fc Steidell 
(|2009l ). Therefore we suggest that more than half of total 
o"DLA are hosted by the fainter LBGs than the current ob- 
serva tional limit, which is consistent with the observation 
(e.g.. lFvnbo. MoUer fc Warren I [l999l ). 

The green shaded region is an approximate range 
of SFR for DLAs derived by the Cfnj cooling technique 
j Wolfe. Prochaska fc GawiseJ |2003| ). More than eighty per 
cent of DLA cross section in our simulations overlaps with 
the green shaded range. There are some observed DLA 
galax;ies with higher SFR of a few tens MQyr~^, which 



are related to LBGs (e.g., iMoUer et all |2002| ). The mas- 
sive haloes with A'/h > 10"'^ M© in our simulation would 
correspond to those DLA host galaxies with relatively high 
SFR. 

The bottom two panels of Figure [11] show the PDF as 
functions of stellar mass and gas metallicity. The median 
stellar mass is ~ 2.4 x 10* M©, and the median metallicity 
is {Z/Zq) = 0.1. Fifty per cent of total ctdla are hosted 
by galaxies with 3.9 x lO'' < K'L/Mq < 1.8 x 10^ and 
0.06 < Z/Zq < —0.2, both centered around the median val- 
ues (i.e. ±25% around median). The distribution in terms of 
metallicity is strongly peaked around (Z/Zq) — 0.1, and has 
an extended tail to lower metallicities. Our mean/median 
metallicity is not s o far from the peak o f observational 
data, Z/Zq ~ 0.1 (jProchaska et al]|2007bl ). However, our 
PDF is somewhat shifted to hig her metallicity than the 
observed DLA's at 2 ^ z ^ 4 (IProchaska et al.1 [2007bl ) 
and 1.5 ^ z < 5 (|Rafelski et al.1 l2012t ). Our simulation 
can successfully reproduce f{Nm), but it seems to over- 
predict DLA metallicity distribution. The metallicity-A^'Hi 
rel ation is not sensitive to the wind models (e.g.. Figure 9 
in iNagamine et al.ll2004lj ). Hence, our current wind model 
does not seem to cause the higher metallicity. The higher 
metallicity may be caused by the star formation model. If 
our star formation model is too efficient at high redshift, 
the gas is quickly enriched by metals, and it leads to higher 
metallicity of DLAs. In addition, the resolution of numer- 
ical simulations may change the metallicity distribution of 
DLAs. In general, higher resolution simulations can resolve 
smaller dwarf galaxies, and many gas clumps in them. If 
the gas in the dwarf galaxies which have little metals can 
contribute to the DLA abundance, the typical metallicity 
of DLAs can be reduced. These topics will be addressed in 
future works. 

Another possibility is the selection bias of DLAs in 
observation due to the dust r eddening (e.g.. [ Ellison et al] 
200ll:lFynbo et al. | l20ld l201ll : [Khare et al.|[2012h . Recently 



Fynbo et al.] (|2010l . l201ll ) found some metal-rich DLAs using 



the detection technique based on Sill and Fell absorption 
lines. These objects were not detecte d in the automatic DLA 
searches in the SDSS using Lya (e.g.. [Prochaska et al.|[2005l ). 
They suggested that dust in such metal-rich DLAs could ob- 
scure the QSOs, and cause the non-detection of DLAs using 
Lya. This would lead to a distribution of DLAs biased to 
low metallicities. Hence, the difference between our simula- 
tions and observat ion may come from this selection bias. 
[Khare et al] (|2012t ) concluded that 10 per cent of the 
DLAs in SDSS DR7 cause significant reddening. It is unclear 
if this reddened population can explain the discrepancy we 
see in the PDF, but it will bring the two distributions closer. 

The dashed histograms are the PDFs in the optically 
thin run. Although f{Nui) significantly changes with the 
UVB model, the typical physical quantities do not change 
very much. As we see in Figure [31 P in the fitting function 
(Equation 5) changes with the UVB model, while the power- 
law slope a does not change largely. In our simulations, the 
physical quantities are roughly proportional to halo mass. 
Therefore, the typical physical quantities do not depend on 
the UVB model very much. 

Figure [H shows the PDF of ctlls. The shape of PDF 
is basically similar to that of ctdla. However, in the case 
of LLSs, the contribution from lower mass (and hence 
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Figure 11. The PDF of DLA cross section in the OTUV run 
at z = 3 as a function of various physical quantities: halo mass, 
SFR, stellar mass, and metallicity, from top-left to lower right, 
respectively. The dashed histograms are PDFs in the Optically- 
Thin run. The vertical dashed and dotted lines are the median 
and moan of PDF, respectively. The green shaded region is the 
halo mass ra nge of DLAs host g alaxies suggested by the cluster- 
ing analysis l lCooke et al.|[2006h . The cyan shaded region is an 
approxima te range o f SFR fo r DLAs derived by the C[ii] cooling 
technique l lWolfe. Procliaska fc Gawise i1 l2003l) . The vertical, blue 
dot-dashe d line in the upper righ t panel is the SFR, of L*— LBGs 
at = 3 ijReddv &: Steideill2009l ). The purple histogram in the 
bottom right p anel shows the PDF of ob served DLA metallicities 
at 2 < z ^ 4 bv lProchaska et al.l l|2007bl '). The blue dashed curve 
is the Gaussian fit to the metallicity of 207 DLAs at 1.5 ^ z ^ 5 
dRafelski et al.ll2012t) . 







(bias) 


(SFR) 


(A/.) 


(Z/Zq) 




(IQio A/q) 




(Moyr-i) 


(10« A/q) 




DLA 


2.39 


1.69 


0.33 


2.43 


0.10 


LLS 


0.96 


1.50 


0.06 


0.65 


0.08 



Table 2. Median values of various physical quantities for 
the PDF of DLA/LLS cross section. The bias is derived from 
the median halo mass by using an ellipsoidal collapse model 
llSheth. Mo fc Tormenll200ll) . 



lower SFR and lower M*) haloes becomes greater than 
for DLAs. The median values are: Mh = 9.6 X 10^ Mq, 
SFR = 0.06 MQyr^\M* = 6.5 x 10^ A-Zq , and Z/Zq =0.08. 
About 80 per cent of total ctlls come from haloes with 
SFR < 1 M0yr~^. Hence the majority of LLSs are hosted 
by the low-mass halos with LBGs that are fainter than the 
current observational limit. Fifty per cent of total ctlls are 
hosted by haloes with 2.5 x 10^ < M^,/Mq < 5.3 x 10** and 
0.05 < Z/Zq < 0.15. We summarize the various median 
values in Table [2] 



Figure 12. Same as Fig. 1111 but for LLS cross section. 



4 SUMMARY 

We examined the physical properties of DLAs and LLSs at 
z = 3 using cosmological SPH simulations. In particular, we 
studied the impact of radiative transfer of UVB and local 
stellar radiation on DLAs and LLSs. Our major findings are 
as follows: 



with a self-shielding density 
3„„-3 ^jjg UVB can re- 



• The OTUV simulation 
threshold of nl{^ — 6 x 10 " cm 
produce the observed f{Nm) very well (Figure com- 
pared to the run with optically-thin approximation (the 
Optically-Thin run). W e have presented this result earlier 
(|Nagamine et aLlbOlol ). and in this paper we further vali- 
dated the adopted value of n[ji^ by the direct RT calcula- 
tions. Similar resu lts have been obtained various authors to 
a varying degre e dPontzen et al" 20081: Tescari et al. 20091: 



Fumagalli et all 12011 
201lf ) 



McQuinn. Oh fc Faucher-Giguerel 



As we showed in Figs. (2] & El the above value of n[Jj^ 
captures the rapid change in the ionization fraction in halo 
gas very well. In the Optically- Thin run, the UVB sinks too 
deeply into the halo gas and ionizes gas too much, resulting 
in the underestimate of f{Nm)- Given that it is still diffi- 
cult to perform RT calculations simultaneously with hydro- 
dynamics in current cosmological simulations, our OTUV 
model provides a useful prescription for any cosmological 
simulations that cannot resolve the self-shielding by molec- 
ular gas on sub-kpc scales. 

• We find that the local stellar radiation does not strongly 
affect the distribution of DLAs/LLSs, and f{Nm) does not 
change very much by the stellar RT (Figure [T|. This is 
because clumpy high-density clouds near young star clus- 
ters effective ly absorb most of t he ionizing photons from 
young stars (jYaiima et aLllioilal ). The effect of stellar RT 
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is stronger on the lower mass halos, and the reduction in 
o"DLA and ctlls were shown in Figs. [S] 

• We compared the DLA/LLS cross sections in simula- 
tions with different UVB models (Figure |4|. In our simula- 
tions, the No-UVB run gives the highest DLA/LLS cross sec- 
tions as expected, and it overpredicts the rate-of-incidence 
and the column density distribution. For the DLAs, the 
OTUV run gives intermediate (Tdla values, and the HMO. 5 
& Optically-Thin run give the lowest (tdla. The similarity 
of the latter two runs tells us that the gas becomes highly 
ionized under the optically thin run, no matter what the ex- 
act UVB amplitude is. For the LLSs, the latter three runs 
(OTUV, Opticall-Thin, and HMO. 5) all give similar ctlls, 
and the self-shielding model does not matter for LLSs, as 
they are in the outskirts of galaxies. 

• We compared the results of DLA cross section from var- 
ious authors and simulations (Figure[S| . There are some gen- 
eral agreement in the sense that most results are bracketed 
by our Opticall-Thin run and the No-UVB run, and that 
both (Tdla and ctlls increases with increasing halo mass. 
However, there is still significant scatter between different 
authors. The origin of the differences could be a combination 
of many things, such as different models of star formation, 
galactic wind, UVB and self-shielding treatment. For exam- 
ple, the stronger wind model can sweep high-density gas 
clouds out of star-forming regions, and decrease DLA cross 
section. The self-shielding model of UVB can increase the 
DLA cross section compared to the optically-thin method, 
and rate of increase can be different according to the thresh- 
old density for self-shielding. We have chosen the appropri- 
ate value of TL^^ based on the good match to the observed 
f{Nm) and through the validation by direct radiative trans- 
fer calculation. 

• The cumulative number of DLAs per unit redshift 
strongly depends on the UVB models (Figure llOf) . The 
OTUV run agrees well with the observed data. On the other 
hand, the No-UVB run overproduces DLAs. The Optically- 
Thin and HMO. 5 runs underpredict the observation. 

• The visual appearance of DLAs and LLSs is quite dif- 
ferent in the OTUV and Optically-Thin runs (Figure [7] & 
[UJ. In the OTUV run, the DLAs have an extended distribu- 
tion around star-forming regions on typically ~ 10 — 30 kpc 
scales, and LLSs are surrounding DLAs on ~ 30 — 60 kpc 
scales. As shown in the middle row of Figure the extent 
of DLAs are typically smaller than the virial radius of the 
halo, and the LLSs are about the same size as the virial ra- 
dius. Given the good agreement between the OTUV run and 
the observations of /{Nni), we consider that these are the 
most realistic distribution of DLAs and LLSs in our simula- 
tion set. 

On the other hand, in the Optically- Thin run, the UVB 
penetrates deeply, ionizes the hydrogen gas effectively, and 
reduces the column density significantly. As a result, most of 
the region around star-forming galaxies turn into LLSs, and 
the DLAs shrink and become compact. This striking effect 
of optically-thin approximation is highlighted by the com- 
parison of columns (c) and (d) in Figure [71 suggesting the 
importance of self-shielding effect. The pres e nt work super- 
sedes the previous work of iNagamine et al] (|2004bl '). which 
reported clumpy distribution of DLA gas in haloes at z — 3. 
With the proper treatment of self-shielding as in the OTUV 



run, the DLAs have more extended distribution around star- 
forming galaxies on 10 — 30 kpc scales. 

We quantified the above features in the PDFs of impact 
parameters between DLAs/LLSs and the nearest galaxies. 
As expected, we find a much wider distribution for the 
OTUV run than for the Optically-Thin run. 

• We examined the PDF of DLA/LLS cross sections as 
functions of halo mass, SFR, stellar mass, and metallicity 
(Figures [TT] and I12p . The median values of distributions 
are summarized in Table [2] About 30 per cent of DLAs 
are hosted by galaxies having SFR =1 — 20 Moyr"^, 
which is the typical SFR range for LBGs. More than half 
of DLAs ((Tdla, to be more precise) are hosted by the 
LBGs that are fainter than the current observational limit. 
About 80 per cent of total ctlls are hosted by haloes with 
SFR < 1 Moyr"^, hence most LLSs are associated with 
low-mass halos with faint LBGs below the current detection 
limit. 
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